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1 Introduction 

Bound states in integrable quantum field theory models in 1+1 dimensions have been 
studied extensively for many years [1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12]. The quantum 
bound states are usually constructed by using either the coordinate Bethe ansatz or the 
algebraic Bethe ansatz. For an integrable non-relativistic Hamiltonian, the coordinate 
Bethe ansatz can yield the exact eigenf unctions in the coordinate representation. If such 
an eigenfunction decays sufficiently fast when any of the particle coordinates tends towards 
infinity (keeping the center of mass coordinate fixed), we call such a localized square- 
integrable eigenfunction a quantum bound state. It is also possible to construct quantum 
bound states using the algebraic Bethe ansatz, by choosing appropriate distributions of 
the spectral parameter in the complex plane [3, 4, 5]. It is usually found that localized 
quantum bound states of various integrable models, including the well known nonlinear 
Schrodinger model (NLS) and the sine-Gordon model, have positive binding energy [1, 2, 



In this paper, we will study the quantum bound states of a generalized NLS model. 
Classical and quantum versions of this type of models have found applications in different 
areas of physics like circularly polarized nonlinear Alfven waves in a plasma [13, 14] and 
quantum properties of solitons in optical fibers [15]. The Hamiltonian of the generalized 
NLS model in its second quantized form is given by [7] 



where the subscripts t and x denote partial derivatives with respect to time and space 
respectively, rj and c are real coupling constants, and we have set the particle mass m = 
1/2. The coupling constant i] is dimensionless whereas c has the dimension of inverse 
length. The field operators ip(x,t), ip^(x,t) obey the equal time commutation relations, 

ty{x,t),il){y,t)\ = W(x,t),^{y,t)] = 0, and [^(x,t),^(y,t)] =h5{x-y). 

It may be observed that for two special cases corresponding to i] = and c = 0, 
the Hamiltonian of the generalized NLS model (1.1) reduces to that of the NLS model 
and the derivative NLS (DNLS) model respectively. The quantum bound states for these 
integrable NLS and DNLS models have been investigated earlier by using both algebraic 
Bethe ansatz and coordinate Bethe ansatz [1, 2, 3, 4, 7, 8, 9, 10, 11, 12]. In this paper, 
our aim is to investigate the ranges of values of the coupling constants for which localized 
quantum iV-body bound states exist for the generalized NLS model. In this context it 
should be noted that, unlike the cases of the NLS and DNLS models, the integrability 
of the generalized NLS model (1.1) has not been established so far. So the method of 
algebraic Bethe ansatz can not be used at present to study this model. 

It may also be mentioned that an analysis of the classical NLS model suggests a 
fascinating application of the spectral analysis for nonlinear operators, which has been 
studied earlier [16]. However, for the generalized NLS model, we shall take the following 
approach. Rather than trying to deal directly with the spectral analysis of the model 
defined in (1.1), we shall project the second quantized Hamiltonian to the bosonic N- 
particle subspace and reduce it to a linear operator. The method of coordinate Bethe 



3, 4, 5, 6]. 




(1.1) 
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ansatz can then be used to find the spectrum of such a linear operator. In Sec. 2, we 
apply this coordinate Bethe ansatz to find out the conditions which the Bethe momenta 
have to satisfy in order that a quantum iV-body bound state should exist. In Sec. 3, 
we analyze these conditions on Bethe momenta and find that such bound states exist for 
all possible values (both positive and negative) of c and within several non-overlapping 
ranges (called bands) of rj. We also apply the idea of Farey sequences in number theory 
to completely determine the ranges of all bands for which iV-body bound states exist 
for a given value of N. In Sec. 4, we show that the bound states appearing within 
each band can have both positive and negative momentum. Moreover, bound states with 
zero momentum can be constructed for a wide range of the coupling constants. We also 
calculate the binding energy of the bound states and find that it can take both positive 
and negative value within each band. In Sec. 5 we make some concluding remarks. 



2 Conditions for quantum iV-body bound states in 
the generalized NLS model 

To apply the coordinate Bethe ansatz, we separate the full bosonic Fock space associated 
with the Hamiltonian (1.1) into disjoint iV-particle subspaces \Sn)- We want to solve the 
eigenvalue equation H\Sn) = E\Sn). The coordinate representation of this equation is 
given by 

H N t n (x 1 ,x 2 , • • • ,X N ) = E T N (X!,X 2 , • • • ,x N ) , (2.1) 
where the iV-particle symmetric wave function r N (xi,x 2 , ■ ■ ■ , xjv) is defined as 

t n (x 1 ,x 2 ,- ■ ■ ,x N ) = -^(0\'4>(x 1 )---'4>(x N )\S N ), (2.2) 

and if 7v, the projection of the second-quantized Hamiltonian if (1.1) on to the N particle 
sector, is given by 

H N = -h 2 £ — 2 + 2h 2 c £ 5{ Xl -x m ) + 2ih 2 r, £ *(*,-*„,) t + (2.3) 

j=l UX j Km Km UX l UX ™ 

It is evident that if at commutes with the total momentum operator in the iV-particle 
sector, which is defined as 

N Q 

p N = -in £ ] -. (2.4) 
j=i ux j 

Note that if jv remains invariant while P/v changes sign if we change the sign of r] and 
transform all the Xi — > —Xi, keeping c unchanged; let us call this the parity transformation. 
Hence it is sufficient to study the model for i] > 0. The eigenf unctions for r] < can then 
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by obtained by changing Xi — > — xf, this leaves the energy invariant but reverses the 
momentum. 

Let us first construct the eigenfunctions of the Hamiltonian (2.3) for the two-particle 
case, without imposing any symmetry property on t 2 (xi,x 2 ) under the exchange of the 
particle coordinates. For the region x\ < x 2 , we may take the eigenfunction to be 

t 2 (x 1 ,x 2 ) = exp {i{kiXi + k 2 x 2 )} , (2.5) 

where k\ and k 2 are two distinct wave numbers. Using Eq. (2.1) for N = 2, we find that 
this two-particle wave function takes the following form in the region X\ > x 2 : 

t 2 (xi,x 2 ) = A(ki,k 2 )exp {i(kiXi + k 2 x 2 )} + B(k 1 , k 2 ) exp {i(k 2 xi + kix 2 )} , (2.6) 

where the 'matching coefficients' A(k±, k 2 ) and B(ki, k 2 ) are given by [7] 

A(k u k 2 ) = '- , B(h,k 2 ) = 1 - A(k 1: k 2 ) . (2.7) 

Ki — K 2 

By using these matching coefficients, we can construct completely symmetric iV-particle 
eigenfunctions for the Hamiltonian (2.3). In the region x\ < x 2 < ■•• < Xn, these 
eigenfunctions are given by [7, 17] 

r N (x 1 ,x 2 ,---,x N ) = ( II ~~77T^ rfi- I Pu(i)*><2),"-rtX)( x i> x 2> ■■■> x n) , (2.8) 

u> \l<m A \ K m-, Kl) J 

where 

Pw(i)M2),"-M*)( x i> x 2>'''> x >r) = ex P {i(k u >(i)X 1 H h k^XN)} ■ (2.9) 

In the expression (2.8), the k n 's are all distinct wave numbers, and XL implies summing 
over all permutations of the integers (1, 2, ....N). The eigenvalues of the momentum (2.4) 
and Hamiltonian (2.3) operators, corresponding to the eigenfunctions t^(xi,x 2 , • • • , xn), 
are given by 

N 

P N t n (xi,x 2 , ■ ■ ■ ,x N ) = kj^j t n (xi,x 2 , ■ ■ ■ ,x N ) , (2.10a) 

3=1 
N 

H N t n (xi,x 2 , ■ ■ ■ ,x N ) = h 2 (Y^ k f) t n(xi,x 2 , ■ ■ ■ ,x N ) . (2.106) 

3=1 

The wave function in (2.8) will represent a localized bound state if it decays when 
any of the relative coordinates measuring the distance between a pair of particles tends 
towards infinity. To obtain the condition for the existence of such a localized bound state, 
let us consider the following wave function in the region x\ < x 2 < • • • < Xn- 

N 

Pi,2,-,jv( x 1 ,x 2 , x N ) = exp (iJ^kjXj) . (2.11) 

3=1 
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As before, the momentum eigenvalue corresponding to this wave function is given by 
hJ2f=i kj. Since this must be a real quantity, we obtain the condition 

Efc = 0, (2.12) 

3=1 

where qj denotes the imaginary part of kj. The probability density corresponding to the 
wave function Pi$,-,n{ %i,%2, • • • ,%n ) in (2-11) can be expressed as 

N-l r 

|Pi,2,...,jv( xi,x 2 ,---,x N )| 2 = exp { 2 E Vr } > (2-13) 

r=l jr'=l 

where the y r 's are the N — 1 relative coordinates: y r = a; r+ i — x r , and we have used Eq. 
(2.12). It is evident that the probability density in (2.13) decays exponentially in the 
limit y r — > oo for one or more values of r, provided that all the following conditions are 
satisfied: 

N-l 

9i <0, gi + g 2 <0, , £ <&<0. (2- 14 ) 

3=1 

Equations (2.12) and (2.14) imply that the wave function Pi i 2,-,at(^i, x 2 , ■ ■ ■ , a; at) (2.11) 
is square-integrable if one holds the centre-of-mass coordinate X = J2i x i/N fixed, and 
integrates over the relative coordinates y r . In the region x\ < x 2 < ■ ■ ■ < x^, the integrals 
over the y r 's all run from to oo, and they are independent of each other. The probability 
density (2.13) is independent of X and due to the conditions (2.14), the integration of 
this probability density over the y r 's gives a finite result. 

Note that the wave function (2.11) is obtained by taking uj as the identity permutation 
in (2.9). However, the full wave function (2.8) also contains terms like (2.9) with uj 
representing all possible nontrivial permutations. The conditions which ensure the decay 
of such a term with a nontrivial permutation will, in general, contradict the conditions 
(2.14). To construct a decaying wave function, therefore, the coefficients of all terms 
Pu)(i),u)(2),-,u)(N)(xi, X2, • • ■ , Xn) with nontrivial permutations must be made to vanish. It 
turns out that it is sufficient to require: 

A(h,k 2 ) = 0, A(k 2 ,k 3 ) = 0, , A(k N _ u k N ) = 0. (2.15) 

Thus the simultaneous validity of the conditions (2.12), (2.14) and (2.15) ensures that 
the full wave function t^(xi,x 2 , ■ ■ ■ ,xn) (2.8) represents a localized bound state. Using 
the conditions (2.12) and (2.15), one can obtain an expression for the complex fc n 's in the 
form 

K = + (2.16) 

2 tan0 

where x is a rea l parameter, and is related to the coupling constant as 

= tan -1 (?7) . (2.17) 

To obtain an unique value of from the above equation, we restrict it to the fundamental 
region — | < 0(^ 0) < |. [Note that rj and <fi have the same sign. Due to the parity 
symmetry mentioned above, we can restrict our attention to the range < < |]. 
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Now, let us verify whether the fc n 's in (2.16) satisfy the conditions (2.14). Summing 
over the imaginary parts of these fc„'s, we can express the conditions (2.14) in the form 

^Sm(Z0) si r/ 7V _An > q for l = x 2> ... N _ 1 / 2 _ lg x 

sin 

Thus, for some given values of c, 0, TV and x, a bound state will exist when all the above 
inequalities are simultaneously satisfied. Surprisingly, the above condition for having 
a quantum bound state for the generalized NLS model is completely independent of 
the coupling constant c. By using Eqs. (2.16) and (2.10a), we obtain the momentum 
eigenvalue of such a bound state to be 

- - tlCN + (2.19) 



2 tan sin 
Hence x is related to the momentum as 

v = h P -™L) Sin< ^ . (2.20) 
X h { 2tan0 ; S in(iV0) 1 ; 

In a similar way, the energy eigenvalue of these bound states is found to be 

ft V sin(2iV0) h 2 c 2 N ft 2 xcsin(jV0)cos0 

sin(20) + 4tan 2 + sin 2 ' 1 } 

For some fixed value of c, and N, the above expression for the energy has a minimum 
for a nonzero value of the momentum given by 

_ iVftc tan(iV0) 
~ 2tan0^ ~ iVtan0 ' ' 

provided that tan(iV0) > 0. The next section of our paper will be devoted to finding the 
ranges of values of where all the inequalities (2.18) are simultaneously satisfied for a 
given value of the particle number N. 



3 Determining the values of where TV-body bound 
states exist 

In this section, we will study the values of where iV-body bound states exist for different 
values of iV for the generalized NLS model. In an earlier analysis, the quantum bound 
states of this model had been found to exist for < < [7]. However from our 
analysis, it will turn out that there are several non-overlapping ranges (bands) of for 
which quantum bound states exist for the generalized NLS model. Even the lowest band 
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for which such bound states exist (i.e., when lies in the range < < j^i) is wider 
than that obtained earlier [7]. Since eqn. (2.18) is independent of the coupling constant 
c, bound states are formed within the above mentioned bands for all possible values of 
c. Consequently, in contrast to the case of NLS model where bound states appear only 
for negative values of c, bound states of the generalized NLS model can appear for both 
positive and negative values of c. 

Incidentally, it is quite remarkable that our results concerning the band structure of 
the generalized NLS model turn out to be very similar to our analysis on the quantum 
bound states of the DNLS model [11, 12]. For the simplest case N = 2, the condition 
(2.18) is satisfied when lies in the range < < f (— f < < 0) for the choice x > 
(x < 0). Thus any nonzero value of within its fundamental region can generate a 2-body 
bound state. 

We will now consider the more interesting case with N > 3. Due to the parity 
symmetry of the Hamiltonian in (2.3), we will henceforth assume that > 0. 

Let us consider a value of of the form 

nn 

(PN,n = -Jjr , (3.1) 

where n is an integer satisfying 1 < n < N/2. 

It can be easily shown that all the inequalities in (2.18) are satisfied for = <pN,n, if 
and only if N and n are relatively prime (with n odd for x > 0, and n even for x < 0) 
[12]. By continuity, it then follows that all the inequalities will hold in a neighborhood of 
0_/v,n extending from a value 4>N,n,- to a value 4>N,n,+, such that 0jv,n- < 4>N,n < <t>N,n,+- 
We will call the region 

0AT,n - < < 4>N,n,+ (3.2) 

as the band B^ n . In this band, there is a bound state with N particles. 

We now have to determine the end points 0jv,n,- and <pN,n,+ of the band Bn,u- The 
inequalities in (2.18) show that the end points are given by of the form 

= y , (3.3) 
where j and I are relatively prime and satisfy the conditions: 

1 < I < N , and j < -. (3.4) 

Thus the end points of the band Bm, u are given by two rational numbers <f)/ir of the form 
j/l which lie closest to (and on either side of) the point 0Ar 5 „/vr = n/N . The solution to 
this problem is well known in number theory and is described by the Farey sequences [18] . 

For a positive integer N, the Farey sequence F N is defined to be the set of all the 
fractions a/b in increasing order such that (i) < a < b < N, and (ii) a and b are 
relatively prime. 
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The Farey sequences for the first few integers are given by 

1 



1 1 

1 1 

12 1 

1 1 



1 3 2 3 1 

1112 3 1 

1 4 3 2 3 4 1 

1 1 1 2 1 3 2 3 4 1 

15435253451 



(3.5) 



To return to our problem, we now see that the points <pN,n in (3-1) (which lie in the 
bands B N;n ) have a one-to-one correspondence with the fractions n/N, which appear on 
the left side of 1/2 within the sequence F N . Due to Eqs. (3.3) and (3.4), the end points 
of the band B^^ are given by 

TTCli 710,2 
<PN,n- = -r- , <PN,n,+ = -7— , (3.6) 

where a,\/h\ and 02/^2 are the fractions lying immediately to the left and right of n/N in 
the Farey sequence F N . These are the two unique fractions which lie closest to (and on 
either side) of n/N. 

In Table 1, we show the ranges of values of for which bound states exist for N = 2 
to 9. 



N 


n 


Range of values 






of 4>/tt 


2 


1 


< 0/tt < 1/2 


3 


1 


< (f)/n < 1/2 


4 


1 


< 0/vr < 1/3 


5 


1 


< 0/vr < 1/4 


5 


2 


1/3 < <p/ir < 1/2 


6 


1 


< 0/vr < 1/5 


7 


1 


< 0/vr < 1/6 


7 


2 


1/4 < 0/tt < 1/3 


7 


3 


2/5 < 0/tt < 1/2 


8 


1 


< 0/tt < 1/7 


8 


3 


1/3 < 0/tt < 2/5 


9 


1 


< 0/tt < 1/8 


9 


2 


1/5 < 0/tt < 1/4 


9 


4 


3/7 < (p/rr < 1/2 



Table 1. The range of values of (p/n for which bound states exist for various values of N. 



8 



4 Momentum and binding energy of a TV-body bound 
state 



In the previous section, we have determined the band structure associated with the gen- 
eralized NLS model by proceeding in the same way as in the case of the DNLS model. 
However it should be noted that, in the case of the DNLS model, one does not get quan- 
tum bound states with zero momentum expect for some discrete values of the coupling 
constant = nn/N [12]. Thus, in general, bound states of the DNLS model do not exist 
in the centre-of-mass (COM) frame. This may be due to the fact that in the Hamiltonian 
of the DNLS model there is no dimensionful coupling constant, and it is invariant under a 
scale transformation. On the other hand, the Hamiltonian of the generalized NLS model 
(1.1) has a dimensionful coupling which breaks the scale invariance. So, we may expect 
to get zero momentum bound states for this model within a wide range of the coupling 
constants. The binding energy of such zero momentum states may then be interpreted 
as their internal energy in the COM frame. With this aim in mind, in this section we 
will calculate the momentum and binding energy for the iV-body bound states described 
above. 

We first look at the momentum of the bound states in a particular band B Nn using 
Eq.(2.19). The form of the end points given in (3.6) shows that sin(iV0) = at only one 
point in the band B^ n) namely, at = <pN,n- in the right part of the band (i.e., from <pN,n 
to 4>N,n,+)i the sign of sin(iV0) is (—1)™. In the left part of the band (i.e., from <f>N, n ,- 
to 4>N,n)i the sign of sin(iV0) is (— 1)™ +1 . Now, the analysis given above showed that x 
has the same sign as (— l) n+1 . Hence the momentum given in (2.19) is greater (less) than 
hcN / (2 tan0) in the left (right) part of a band, and is equal to UcN / (2 tan0) at = 4>N,n- 

In Figs, la and lb, we have plotted the allowed values of momentum (P) (2.19) as a 
function of (p/n for N — 5. In Fig. la, (taking he = 1) the minimum (maximum) value 
of P in the left (right) part of a band is 5/2 tan 0. Similarly, in Fig. lb (taking he = — 1) 
the minimum (maximum) value of P in the left (right) part of a band is —5/2 tan 0. 

Next, we look at the energy using Eq. (2.21). To calculate the binding energy, we 
consider a reference state in which the momentum P of the iV-body bound state given in 
(2.19) is equally distributed amongst N single-particle scattering states. The real wave 
number associated with each of these single-particle states is denoted by k$. From Eqs. 
(2.10a) and (2.19), we obtain 

k - c Xsin(jV0) 
h ~ 2W + iVsin0 ' 

Using Eq. (2.10b), we can calculate the total energy for the N single-particle scattering 
state as 

E , = n-Nkl = »V™ a W> + + fiVri-WW, (4 . 2) 

TV sin 2 4tan 2 shr0 

Subtracting E in (2.21) from E s in (4.2), we obtain the binding energy of the iV-body 
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bound state as 



E M ,N) S E.-E = »V*W>jg(M _ . (4.3) 

sin0 L Nsm<p cos0 J 

It may be noted that the above expression for the binding energy remains invariant under 
the transformation — + — 0. Now, it can be easily shown that this expression for the 
binding energy is positive in the left part of each band, negative in the right part, and 
zero at the point = <pN,n [12]. 

We thus see that for > 0, (P — /iciV/2tan0) and the binding energy are both 
positive in the left part of a band, and they are both negative in the right part of a band. 
Consequently, for the case > and c > (c < 0), bound states with zero momentum 
can be constructed within the right (left) side of each band. Using eqn. (2.18) it is easy 
to see that, the momentum of a quantum bound state vanishes for 

cN cos0 , . 

X = ~~2~sin(iV0)' ( • ' 

Substituting this value of x to eqn. (4.3), we obtain the binding energy for these zero 
momentum bound states as 

^p-n, , h 2 c 2 N r iVtan0 1 . . 

E£-°(6,N) = - 5-r 1 • 4.5 

B vr ' ; 4tan 2 L tan(iV0)J v ; 

This can be interpreted as the internal energy of a quantum bound state in its COM 
frame. 

In Fig. 2a and 2b, we have plotted the above expression for the binding energy (4.5) 
as a function of (f)/ir for iV = 5 (he = 1 in Fig. 2a and he = — 1 in Fig. 2b). 

Finally, we want to discuss a subtle connection between quantum bound states of the 
generalized NLS model and the ordinary NLS model. It has been already mentioned that, 
in the limit (f) — > 0, the ordinary NLS model can be obtained from the generalized NLS 
model. It is well known that quantum bound states of NLS model exist only for c < 
with the iV momenta being given by 

kn = ^ + i\ (N+l-2n) , (4.6) 

and the binding energy for such state taking the form 

E B (N) = h 2 c 2 N(N 2 - 1)/12. (4.7) 

On the other hand, in this article we have seen that, for any small value of the parameter 
lying within the lowest band, quantum bound states of the generalized NLS model exist 
for both c > and c < 0. So it is natural to ask why bound states corresponding to c > 
disappear in the limit 0^0. To answer this question, we first try to take the — » limit 
of the complex momenta k n (2.16) associated with the generalized NLS model and check 
whether they reduce to the k n given in eqn. (4.6). However it is evident that, all k n in 
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eqn. (2.16) as well as the corresponding total momentum P (2.19) will diverge if we take 
the limit — > by keeping x fixed. To bypass this problem, we have to take the limit 
— > in a way such that the total momentum P (2.19) remains constant. For the sake of 
convenience we choose P = 0, for which the value of x is given in eqn. (4.4). Substituting 
this value of x i n ec l n - (2-16) and expanding the right hand side of this equation in a 
power series of 0, we easily obtain 

k n = i C -{N + l-2n) + 0(0), (4.8) 

which reproduces the value of k n in eqn. (4.6) for P = 0. Similarly, by taking the limit 
— > of the binding energy (4.5) corresponding to the zero momentum bound states of 
the generalized NLS model, we can reproduce E B {N) in eqn. (4.7). In the case of the 
generalized NLS model with c > 0, we have already observed that zero momentum bound 
states exist only in the right side of the each band. For taking the limit — » 0, however, 
we have to move to the left side of the lowest band. Consequently, it is impossible to take 
the limit — > within the lowest band after fixing P = (same conclusion can be drawn 
for any other fixed value of P). On the other hand, for c < 0, zero momentum bound 
states exist in the left side of each band and one can easily take the limit after 
fixing P = 0. This fact clearly explains why quantum bound states remain in the limit 
0^0 only if c < 0. 



5 Conclusion 

By applying the coordinate Bethe ansatz, we have investigated the range of the coupling 
constants (c,0) for which localized quantum iV-body bound states exist in a generalized 
NLS model. It is found that such bound states exist for all possible values (both positive 
and negative) of c and within several non-overlapping ranges (called bands) of 0. Using the 
ideas of Farey sequences appearing in number theory, we have given explicit expressions 
for all the allowed bands of for which iV-body bound states exist. 

We have also calculated the momentum and binding energy for the bound states within 
all bands in the region > 0. It is found that bound states with positive (negative) values 
of (P — hcN / 2 tan. (f>) appear in the left (right) part of each band, and these states have 
positive (negative) binding energy. Consequently, bound states with zero momentum can 
be constructed within a wide range of the coupling constants. This happens in contrast 
to the case of the DNLS model, where zero momentum bound states appear for only 
some discrete values of 0. The binding energy of zero momentum states has also been 
calculated and interpreted as the internal energy in the centre-of-mass frame. We have 
also remarked on some subtleties of taking the limit — > for the bound states of the 
generalized NLS model, and explained why this process does not lead to any bound state 
of NLS model for c> 0. 
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It is currently not known if the generalized NLS model discussed here is integrable. 
If the model turns out to be non-integrable, then the bound states which we have found 
using the Bethe ansatz may not be complete, and there could be some other bound states. 
This may be an interesting problem for future studies. 
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Fig. 1. The filled in regions indicate the values of the momentum as a function of 4>/ti for 
which a 5-body bound state exists. In figures (a) and (b), hc= 1 and —1 respectively. 
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Fig. 2. The binding energy of the zero momentum 5-body bound state as a function of 
(p/n. In figures (a) and (b), Hc—1 and —1 respectively. 
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